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Strong spin-orbit coupling (SOC) effects of heavy d-orbital elements have long been neglected in 
describing the ground states of their compounds thereby overlooking a variety of fascinating and 
yet unexplored magnetic and electronic states, until recently. The spin-orbit entangled electrons in 
such compounds can get stabilized into unusual spin-orbit multiplet J-states which warrants severe 
investigations. Here we show using detailed magnetic and thermodynamic studies and theoretical 
calculations the ground state of Ba 3 ZnIr 209 , a 6H hexagonal perovskite is a close realisation of the 
elusive J — 0 state. However, we find that local Ir moments are spontaneously generated due to 
the comparable energy scales of the singlet-triplet splitting driven by SOC and the superexchange 
interaction mediated by strong intra-dimer hopping. While the Ir ions within the structural Ir209 
dimer prefers to form a spin-orbit singlet state (SOS) with no resultant moment, substantial in¬ 
terdimer exchange interactions from a frustrated lattice ensure quantum fluctuations till the lowest 
measured temperatures and stabilize a spin-orbital liquid phase. 


The electronic and magnetic properties of 5d transi¬ 
tion metal (TM) compounds often exhibit unusual prop¬ 
erties due to the presence of strong spin-orbit coupling 
(SOC), which is comparable to their on-site Coulomb 
([/) and crystal field (Acfe) interactions [T]. While 
quantum numbers mi (orbital) and mg (spin), for each 
electron can still be considered good and the spin-orbit 
coupled many-body multiplets can be denoted by spec¬ 
troscopic term symbols in the weak spin-orbit (L-S) 
coupling limit, the representation gets strongly modified 
for heavy multielectron elements having strong SOC. In 
such j-j coupling regime, total Mj {X,mj) becomes the 
only valid quantum number and the multiplets and their 
degeneracies are all solely determined by the total an¬ 
gular momentum J. 

The electronic and magnetic response of a system in 
the j-j coupling limit are not yet well-understood and 


therefore, in the recent time, has generated significant 
curiosity. For example, the unconventional insulating 
state of the layered tetravalent iridates (Ir"‘+; 5d®), des¬ 
ignated as novel jeff = | Mott insulators, have at¬ 
tracted much interest HHg, as metallic behavior is usu¬ 
ally expected for partially filled broad t 2 g bands in iri¬ 
dates. The insulating behavior is explained within single 
particle theories by assuming splitting of t 2 g bands into 
a fully filled Je// =3/2 quartet bands and half filled 
narrow Jeff = 1/2 doublet bands due to finite SOC, 
which further splits into a fully occupied lower Hub¬ 
bard and empty upper Hubbard band in the presence of 
relatively small Hubbard U. 

Another intriguing case would be a pentavalent Irid¬ 
ium compound (Ir®+; 5d^), where all the spin-orbit en¬ 
tangled electrons will be confined to the singlet J = 0 
{Mj = 0) state, without any resultant moment or mag- 
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netic response. This is shown in Fig. 1(a) where the 
energy levels of a low spin 5d Ir^+ ion calculated as a 
function of spin-orbit coupling parameter X' [1] evolves 
from the L-S coupling to the j-j coupling regime [6] and 
a J = 0 ground state is realized under strong SOC (see 
Supplementary Information (SI) for details). Recently 
there has been a theoretical suggestion of an excitonic 
magnetism (transition from the J = 0 to J = 1 state) 
in such systems [7], due to comparable energy scales 
of the singlet-triplet splitting determined by SOC and 
the super-exchange interaction promoted by hopping. 
Therefore, if d'^ Ir compounds indeed fall into the cate¬ 
gory of the elusive j-j coupled systems, experimentally 
a zero or vanishingly small magnetic moment should 
be observed for the otherwise nominal S' = 1 Ir^’*' cen¬ 
ter and it would be interesting to know the nature of 
their magnetic response(s), if any. However, stabiliza¬ 
tion of such an unusual electronic state may be subtle 
as this J = 0 ground state at moderate A' could be ex¬ 
ceptionally fragile to minute external perturbation and 
enhanced magnetic responses may appear with rather 
small variations in magnetic fields or IrOg octahedral 
distortions giving rise to noncubic crystal field mi- 

Possibly for this reason, there has been no experimen¬ 
tal realization of this unusual state till date, other than 
a other than a speculation in the case of NalrOs [9], 
while another possible candidate Sr 2 YIr 06 double per- 
ovskite acquired long range magnetic order with com¬ 
paratively larger paramagnetic moment of 0.91 /r^/Ir 
(ideal moment = 2.83 ^b) because of noncubic IrOe 
crystal field [5] . 

In contrast to previous attempts, we have explored 
the complex electronic and magnetic interplay of Ir®"*" 
within a 6H hexagonal compound, namely Ba 3 ZnIr 209 
(Space groupiPOa/mmc; See Fig. 1(b)), where instead 
of the octahedral corner sharing of the IrOe units, the 
octahedra share face and form Ir 2 09 dimers. Our exper¬ 
imental study indeed reveal that Ba 3 ZnIr 209 probably 
comes closest to the coveted J = 0 ground state with a 
low magnetic moment of individual Ir atoms in the para¬ 
magnetic regime, indicating presence of rather strong 
A', but not strong enough to completely arrest the hop¬ 
ping assisted excitonic magnetism at low temperatures. 
Interestingly, first principles density functional calcula¬ 
tions indicate that the adjacent Ir octahedra in a dimer 
prefers to be in a spin-orbital singlet (SOS) state which 
should have eliminated any net magnetization from each 
Ir 209 dimeric unit provided this intradimer singlet in¬ 
teraction is strong enough. Instead, sufficient interdimer 
hoppings are present in order to induce quantum fluc¬ 
tuations in the SOS state and due to frustration aris¬ 
ing from the triangular lattice structure in the a6-plane 
(lower panel of Fig. 1(b)) [7] possibility of any long 
range order is prohibited. Finally, detailed muon spin- 
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FIG. 1. (Color Onlme)(a) The energy states of low spin 5d 
tig atom calculated as a function of spin-orbit coupling 
parameter, (b) The crystal structure of Ba 3 ZnIr 2 09 (upper 
panel) and triangular lattice formed by Ir ions in ab plane 
(lower panel), (c) Experimental (black) and refined (red) 
neutron powder diffraction patterns at room temperature. 


rotation (^SR) and heat capacity measurements reveal 
the presence of strong quantum fluctuations at low tem¬ 
peratures, and absence of any long range magnetic order 
down to 100 mK suggesting Ba 3 ZnIr 209 to be a rare ex¬ 
ample of a quantum spin-orbital liquid (QSOL). 

Polycrystalline Ba 3 ZnIr 209 (BZIO) was synthesized 
by standard solid state reaction using stoichiometric 
amounts of BaC 03 , ZnO and Ir-metal as starting mate¬ 
rials m- The sample purity was checked and refined by 
powder x-ray diffraction measured in the Indian beam¬ 
line (BL-18B) at Photon Factory, KEK, Japan. Neu¬ 
tron powder diffraction (NPD) patterns, recorded in Na¬ 
tional Facility for Neutron Beam Research (NFNBR), 
Dhruva reactor, Mumbai (India) and Paul Scherrer In¬ 
stitute (PSI), Switzerland, both of comparable reso¬ 
lution, were refined by Rietveld method using FULL- 
PROF m- The x-ray photoelectron spectroscopic 
(XPS) measurements were carried out in an Omicron 
electron spectrometer, equipped with EA125 analyzer. 
/iSR experiments were mostly performed with the EMU 
spectrometer at the ISIS large scale facility both in a 
helium flow cryostat and a dilution fridge. A few addi¬ 
tional runs were measured with the DOLLY spectrome¬ 
ter at the PSI continuous source facility where the fast 
early time relaxation at low temperatures can be best 
resolved. The dc magnetic measurements and heat ca¬ 
pacity measurements were carried out using a Quantum 
Design SQUID magnetometer and a Quantum Design 
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PPMS (physical property measurement system), respec¬ 
tively. X-ray absorption spectra at the Zn-itT and Ir 
La-edges were collected at the Elettra (Trieste, Italy) 
ll.lR-EXAFS beamline in standard transmission geom¬ 
etry at room temperature. All the electronic structure 
calculations were carried out using density functional 
theory (DET) within generalised gradient approxima¬ 
tion (GGA) using Hubbard U and SOG as implemented 
in Vienna ab-initio simulation package (VASP) [HIS]. 
This method uses a plane-wave basis set along with pro¬ 
jector augmented waves (PAW) in the ionic core region 

[3 IS]. 

Powder XRD data collected from polycrystalline 
BZIO (not shown) confirm complete phase purity. The 
NPD pattern (with A = 1.48 A; Fig. 1(c)) recorded at 
room temperature has been refined with space group 
PQ^/mmc which revealed < 5% Zn/lr site-disorder, i.e. 
face shared octahedral positions (4/ sites) are almost 
exclusively occupied by Ir forming lr 209 dimer, while 
Zn ions solely occupy the isolated octahedral sites (2a). 
The structural parameters obtained from the refinement 
of the 1.5K data have been listed in Supplementary Ta¬ 
ble I of SI. However, the local cationic disorder and 
octahedral distortions in these compounds may pose a 
different picture |T0], and therefore XAFS experiments 
were carried out. The parameters, obtained from Zn K- 
and Ir L 3 -XAFS data analysis are summarized in Table 
I. The XAFS data analysis confirms negligible chemical 
disorder with Ir (Zn) ions mainly located at the 4/ (2a) 
sites with local interatomic distances being consistent 
with literature. The average Ir-0 distance is compati¬ 
ble with the average crystallographic Ir-0 distance (1.98 
A). It must be noted that the is comparable with 
(and even lower than) the Zn-0 one, which suggests that 
ZnOe and IrOe octahedra have similar distortions. 

Next, the temperature dependence of the electrical re¬ 
sistivity (p) of BZIO, measured by standard four probe 
method, is shown in the inset to Fig. 2(a). The insulat¬ 
ing behavior can be modeled by variable range hopping 
mechanism in 2-dimension as shown in Fig. 2(a). The 
insulating behavior has been further tested by measure¬ 
ment of valence band photoemission experiment, where 
absence of any density of states at the Fermi level (see 
Fig.2(b)) is confirmed. The XPS spectra for the Ir 4/ 
core level (Fig. 2(c)) can be fitted by spin-orbit split 
doublet with separation of 3.04 eV. The energy posi¬ 
tions of 4 / 5/2 and 4 / 3/2 doublet confirm -|-5 oxidation 
state of Ir [11] . 

This intriguing insulating behavior of Ir®"*" having four 
d-electrons in low spin configuration has been investi¬ 
gated within single particle mean-field framework by 
ab-initio density functional theory calculations by sys¬ 
tematically including Hubbard U and SOG. The non¬ 
spin polarized density of states (DOS) of BZIO within 



FIG. 2. (Color Online) (a) The temperature dependence of 
resistivity (inset) can be described by variable range hop¬ 
ping in two dimension (main panel). XPS spectra for (b) 
valence band and (c)Ir 4/ core level, recorded with Mg Ka 
radiation. Total density of states obtained within (d) GGA 
and (e) GGA-I-SOC-I-U. Mechanism of (f) antiferromagnetic 
interaction in presence of SOG and (g) ferromagnetic inter¬ 
action in absence of SOG. 


GGA is shown in Fig. 2(d). As expected, the system 
is metallic with the Fermi level located at the Og man¬ 
ifold, comprising of twelve bands from four Ir atoms in 
the unit cell. The octahedral crystal field is strong, re¬ 
sulting in quite large (^3.5 eV) t 2 g - Cg crystal field 
splitting. The strong intradimer hoppings further lifts 
the degeneracy of the t 2 g states. The DOS is however 
remarkably different upon inclusion of SOG and a Hub¬ 
bard U {U-J= 2.5 eV) which makes the system insu¬ 
lating with a gap 10 meV (see inset of Fig 2(e)). In 
the presence of SOG the t 2 g states (six states including 
spin degeneracy) break into a spin-orbit entangled set 
of states namely two-fold degenerate T 7 and a four-fold 
degenerate Tg states. Since four electrons are available 
per Ir, all the Tg states are occupied and the system is 
insulating as seen experimentally. 

While analyzing the magnetic properties of 
Ba 3 ZnIr 209 in the framework of GGA-flU-SOC 
we find that local moments are spontaneously gener¬ 
ated in the magnetic phase due to the hybridization 
between the occupied Tg and empty Ty states. In 
the following we shall argue that the superexchange 
between a pair of Ir atoms within a dimer will promote 
a spin-orbital singlet state. The four Tg orbitals are 
designated as |la), |l/3), |2a), |2/3) while the two r 7 
orbitals as |3a), |3/3) as illustrated in Fig. 2(f). We 
regard these states as pseudo-spin states, i.e. the 
state in which one electron occupying \ia) (|i/3)) will 
be considered as up (down) pseudo spin states. It is 
interesting to note that unlike t 2 g orbitals the hopping 
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between the pseudo spin states are not always pseudo 
spin conserved m- In particular, an electron in a 
pseudo spin state |1q;) (|1/3)) can hop only to pseudo 
spin states |j/3) (|ja)) where j ^ 1 and all other 
hoppings are pseudo spin conserved m- 

Such hopping leads to antiferromagnetic pseudo spin 
interaction and thereby favors SOS state as shown 
schematically in Fig. 2(f). For an antiferromagnetic 
configuration of pseudo spins in a dimer the Fg elec¬ 
tron in the state lla) and |l/3) can virtually hop to |3/3) 
and 1 3a) respectively as this hopping is not pseudo-spin 
conserved (indicated by in Fig. 2(f). Once these 
electrons are promoted to the r 7 state, they can then 
gain energy by hopping to the other F 7 level as this hop¬ 
ping is now pseudo-spin conserved (indicated by —- in 
Fig. 2(f)). The latter hopping is however not possible 
if the pseudo spins residing in the dimer were paral¬ 
lel. The same dimer however favors ferromagnetic in¬ 
teraction in the absence of spin-orbit interaction as il¬ 
lustrated in Fig. 2(g) as the hopping between the t 2 g 
orbitals are spin conserved. In order to check the above 
mentioned scenarios we have carried out spin polarized 
GGA-I-U and GGA-I-U-I-SOC calculations for several 
magnetic configurations (See SI). The GGA-I-U-I-SOC 
calculations reveal a rather large spin moment (0.65- 
0.75 /iB) per Ir atom. In addition, the large value of 
the orbital moment (0.26 /rB) confirms the importance 
of spin-orbit coupling. Also, the magnetic configuration 
with antiferromagnetic intradimer coupling is found to 
be energetically favorable only in the presence of SOC, 
emphasizing its importance in the realization of SOS 
states in BZIO. Although, an estimate of symmetric ex¬ 
change interactions suggest appreciable intradimer ex¬ 
change (Ji = -14.6 meV), there is substantial interdimer 
exchange too (J 3 = -1.5 meV with 6 neighbors in the 
a — b plane) and (J 2 = -1.6 meV with 3 neighbors per¬ 
pendicular to the a — b plane). Therefore, the intradimer 
interactions are not strong enough compared to the in¬ 
terdimer interactions leading to a frustrated behavior 
instead of a SOS-like long-range state. 

In order to check whether the SOC driven SOSs in 
BZIO eventually lead to long range magnetic order, we 
have carried out NPD experiments at lower temperature 
too. The high resolution NPD pattern recorded at 2 K 
has been shown in Fig. 3(a). All peaks of the observed 
pattern have been refined with PQ^/mmc space group 
without any significant magnetic contribution, indicat¬ 
ing absence of any long range magnetic order in the 
system at least down to 2 K. These NPD results also 
rule out the possibility of any structural transition be¬ 
low 300 K. However, in presence of strong SOC and 
frustration, a clear long range magnetic order might be 
missingabsent while local magnetic interactions are still 
significant. Therefore, in order to unravel the nature of 
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FIG. 3. (Color Online) (a) High resolution neutron powder 
diffraction recorded at 2 K. (b) Time evolution of the muon 
polarization in Ba 3 ZnIr 2 09 in zero external field at different 
temperatures. Lines are fits of the data described in the text, 
(c) Fluctuation rate of the internal fields versus temperature 
as extracted from the fits of the polarizations (see text). 
Inset: At low temperature the fast initial relaxation of the 
muon polarization as measured on the DOLLY spectrometer 
at the continuous muon source of the Paul Scherrer Institut. 


the ground state stabilized in this compound we per¬ 
formed muon spin relaxation experiments in zero field, 
a technique perfectly suited to detect weak and/or par¬ 
tial freezing of local magnetic moments. The evolution 
of the muon polarization in the sample is shown in Fig. 
3(b). At high temperature, the slow Gaussian like re¬ 
laxation arises from nuclear dipoles which are static in 
the /rSR time window and give a T independent con¬ 
tribution. Upon lowering temperature, the relaxation 
rate increases gradually from about 100 K as a result 
of the slowing down of the electronic fluctuations and 
strikingly levels off below about 2 K as demonstrated by 
the nearly perfect overlap of the 2 K and 0.1 K polar¬ 
ization curves. Despite this obvious slowing down of the 
dynamics, we did not observe any spontaneous oscilla¬ 
tions or emergence of a “l/3rd” tail that would signal a 
magnetic transition to a frozen state. The robustness of 
the relaxation to the application of a longitudinal field 
at base temperature also supports its dynamical origin 
(see SI). Therefore, it can be concluded that the system 
does not order down to 0.1 K at least. 

In the low temperature regime, below about 35K, the 
signal appears clearly to be composed of at least 3 com¬ 
ponents; one fast relaxing, one with a moderate relax¬ 
ation rate and one hardly relaxing. These three compo¬ 
nents likely reflect three different muon stopping sites in 
the sample. Therefore we fitted the data to the model. 
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P{t) = ffastG{t,AHi,iy) + f,io^G{t,AH^,v) 

+ (1 - ffast - fslou,)KT{t, AH) ^ ’ 

where KT is the Kubo-Toyabe function m account¬ 
ing for the static nuclear field distribution of width AH, 
ffast and fsiow are the fractions of the two more relax¬ 
ing components and G(t, AH, v) is the dynamical relax¬ 
ation function introduced in Ref. 19 

to account for the effect of a fluctuating field distri¬ 
bution of width AH at a rate v. The fits with respect 
to this model are shown in Fig. 3(b) where only the 
fluctuation rate v was allowed to vary with tempera¬ 
ture and the extracted values are given in Fig. 3(c). 
The other T independent parameters were refined to 
AH = QA{2)G,AHi = 147(5)G,Ai/2 = 39(3)G, 

ffast = fsiow = 31(2)%. At the lowest measured 

temperatures, when the fluctuations are the slowest we 
note that for both relaxing components the inequality 
ly > jAH where 7 = 27r x 135.5 Mrad/s/T is the muon 
gyromagnetic ratio, still holds since yAiJi = 12.5/rs“^ 
and 'fAH 2 = 3.4p,s“^. This validates the use of dynam¬ 
ical relaxation functions, in line with the absence of the 
usual evidence for a frozen ground state. 

The above prediction of a special QSOL state in BZIO 
is further confirmed by the analysis of heat capacity (G). 
The temperature variation of C (Fig. 4(a)) does not 
show any peak or anomaly, consistent with the absence 
of any structural and/or long range magnetic transi¬ 
tion. In order to extract the magnetic heat capacity 
Cm by subtracting the lattice contribution, G was also 
measured for an isostructural nonmagnetic compound 
Ba 3 ZnSb 209 . The difference of molecular weight be¬ 
tween these two compounds for subtracting the lattice 
part has been taken into account following the scal¬ 
ing procedure, developed by Bouvier et. al. [SU]. Cm 
obtained from such analysis has been plotted in Fig. 
4(b), which shows only a broad peak around ^15 K 
followed by a decay around 30 K without any fur¬ 
ther anomaly. A fit of the magnetic heat capacity 
{Cm=lT -\- /3T^) from 3K to 9K gives a significant T- 
linear contribution of 7 = 25.9 mJ/mol-AT^ unusual for 
insulating systems. Also, the finite 7 value does not 
change with application of external magnetic field as 
high as 9 Tesla, suggesting its origin to be gapless ex¬ 
citations from spinon density of QSOL state or minor 
lattice oxygen vacancies and not any paramagnetic im¬ 
purity [2TVl25j . This mixed temperature dependence is 
observed in other spin liquid systems, e.g. Ba 3 CuSb 209 
{S = 1/2) [55], 3C phase of Ba 3 NiSb 209 {S = 1) [57], 
NiGa2S4 (5 = 1) p, Ba3lrTi209 (Jeff = 1/2) [H] 
etc. The magnetic entropy obtained by integrating 
Cm/T with T has been shown in Fig. 4(c), which is 
much smaller (~7%) compared to the value for spin only 



FIG. 4. (Color Online) Temperature dependence of (a) 
specific heat (G), (b) magnetic specific heat (Cm) after sub¬ 
tracting the lattice contribution, (c) Magnetic entropy as 
a function of temperature, (d) FC-ZFC magnetization mea¬ 
sured at 5000 Oe magnetic field fitted using Curie-Weiss law. 
l/(y-yo) vs. T is also plotted showing Ow- 


magnetic entropy, Rln(25-|-1) with 5=1 and unlike 
other spin liquid candidates containing 3d TM, where 
the entropy change is more than 30-50% of the total 
magnetic entropy [5511571 [SO]- The value is in fact 
close to 0 as expected for a spin-orbit coupled mag¬ 
netic entropy, Rln(2J-|-l) with J = 0. Therefore, the 
heat capacity study does confirm presence of QSOL-like 
dynamical ground state in this system, visible below 
30 K. It also indicates towards the fact that the mag¬ 
netic ground state and the moment in BZIO are only 
‘byproducts’ of certain spin excitations from otherwise 
nonmagnetic, singlet Ir®+ J = 0 state, since the small 
magnetic entropy is likely released due to such weak 
magnetic crossover [8]. 

Finally, we show the temperature dependence of field- 
cooled (FC) and zero-field-cooled (ZFC) magnetization 
measured with 5000 Oe field in Fig. 4(d). Clearly there 
is no indication of a magnetic order in this case, as ex¬ 
pected. These x vs. T curves are fit with a function, 
Xo + Cwl{T - Ow) where xo is the temperature inde¬ 
pendent paramagnetic susceptibility, Gw and Ow rep¬ 
resent Curie constant and Curie-Weiss temperature, re¬ 
spectively, and the result is shown as a green straight 
line, laid on the l/(x-Xo) vs. T data. The fitting, espe¬ 
cially at the higher temperature region, extracts a Ow 
value ~-30 K and a moment of ~0.2 gts/li. Clearly, 
the magnetic frustration parameter, defined as the ra¬ 
tio of Ow = -30 K to the upper limit of the magnetic 
ordering temperature that we could probe, going upto 
at least 300 proves that the system indeed behaves like 
a quantum spin-orbital liquid. 









6 


In conclusion, our detailed analysis reveal BZIO is 
a J = 0 Insulator, where the origin of the small but 
finite nonzero moment is due to super-exchange induced 
spin excitations [7]. These Ir moments within the Ir 209 
dimer interact antiferromagnetically resulting in spin- 
orbital singlets. The inter-dimer hoppings in the ab- 
plane induce quantum fluctuations in these singlets and 
a long range magnetic order is arrested due to geometric 
frustration, arising from the triangular lattice structure 
in the a&-plane resulting in a quantum spin-orbital liquid 
state, confirmed both by ^SR and heat capacity studies. 
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1 . 


Derivation of the energy level 


diagram for configuration in L 


S and j — j coupling scheme 


To understand the energy levels of a single state, let us consider the following three-orbital model 
Hamiltonian [1] : 

Htot=HT* + HlO ( 1 ) 

Where, and are respectively the Hamiltonian for the Coulomb interaction and the spin-orbit interaction 
on d orbitals. The Coulomb interaction can be written as. 


Tjint 


Ud ^ + 

2 ^ 1 , 2,3 



i,i' = l,2,3 cy 


Uk 

2 


^ ^ ^ ^ a' 

o-^cr' i,i' = l,2,3 


y + 

i,i' = l,2,3 


( 2 ) 

where, Ud, and Jd are respectively intra-Coulomb interaction, inter-Coulomb interaction and Hund’s rule 
coupling. These Coulomb interactions have the relation Ud = U'^ + 2Jd. dio'(dJo') annihilation (creation) 

operator of the i-th orbital (i = 1,2,3) with a spin a and = dj 


The explicit form of spin-orbit interaction is given by. 


H: 


so 


iX' 

“ y ^' 


E' 


r” ,d' 
cra'^lo 


Imn 


(3) 


where. A' is the magnitude of spin-orbit interaction between orbital {h) and spin (s^) angular momenta of the i-th 
electron and eimn is the Levi-Civita symbol. 

For d'^ filling, the total number of possible configurations are ^C 4 = 15. The eigenvalues and eigenstates are 
calculated by diagonalizing the Hamiltonian Htot in the above mentioned basis set. The variation of the energy 
eigenvalues with A' for U'^ = 1.0 eV and Jd = 0.5 eV is shown in Fig 1(a) of the paper. 

2. Structural parameters after Rietveld refinement of NPD spectrum of Ba-jZnlroOq at 
1.5K.(Table. I) 


3. Details of the Density Functional Theory Calculations 


All calculations reported in this work are carried out using a plane wave based method as implemented in Viena 
ab-initio simulation package (VASP) [HIS]. Exchange and correlation effects are treated using local (spin) density 
approximation (LSDA) with generalized gradient correction (GCA) of Perdew-Burke-Ernzerhof including Hubbard 
U [4] and spin-orbit coupling (SOC). A nominal value of [/ — J = 2.5 eV has been taken unless otherwise stated to 
analyze the effect of correlation on the electronic and magnetic properties of this system. We have used projector 
augmented wave potentials to model the electron-ion interaction. The kinetic energy cut off of the plane wave 
basis was chosen to be 500 eV and 10x10x6 k-mesh has been used for Brillouin-Zone integration. Symmetry has 
been switched off in order to minimize possible numerical errors. 


4. Pseudo spin states 
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The four Tg orbitals are given by, 


and the two r 7 orbitals are, 


|lct) — +*|dzxt))i 


11/3) 

12 a) 


= ~^J^{‘^\dxy\) — \dyz\.) ~ *Mzx^)), 


|2/3) — ^^{2\dxy\) + \dyz\) *Mzxt))l 


|3a) — ^^(\dxy'[) + \dyzi) + i\dzx\))-, 
|3/3) = -^{\dxy{) — Myzt) +*Mzxt)) 


5. Results of spin polarized calculations 


Four different magnetic configurations (see Fig.[^ namely FM (both the intra and inter dimer couplings are ferro¬ 
magnetic), AFMl (both the intra and inter dimer couplings are antiferromagnetic), AFM2 (intra dimer coupling is 
ferromagnetic and inter dimer coupling is antiferromagnetic), and AFM3 (intra dimer coupling is antiferromagnetic 
and inter dimer coupling is ferromagnetic) have been simulated to find out the lowest energy state. Our calcu¬ 
lations as summarized in Table |ll] reveal that AFM2 state has the lowest energy within GGA-l-17 approximation. 
Table 1^ also display the energies of various magnetic states along with spin and orbital moments, obtained within 
GGA-I-C7-I-SOC approach. We find that the SOC changes the lowest energy state from AFM2 to AFMl as discussed 
in the text. 



FIG. 1. The four possible magnetic states. 
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FIG. 2. Asymmetries of the muon decay at 0.1 K measured with various longitudinal fields. 


6. Effect of applied magnetic field on muon spin-rotation (^SR) results 


Fig. 2 shows the evolution of the relaxation when an external field is applied along the initial muon polarization 
at the lowest measured temperature 0.1 K. In case of a static field distribution of width AH^, a strong applied 
field Hl 5AHi should overcome the effect of the internal field and almost completely suppress the relaxation of the 
muon polarization. Clearly the relaxation in the case of Ba 3 ZnIr 2 09 is still strong when a 5000e field is applied 
which ensures the dynamical nature of its ground state like a QSOL state. 
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TABLE I. (Main text)Main structural parameters obtained from the Multi-shell htting of the Zn K-edge and Ir L 3 edge 
XAFS spectra. The crystallographic distances, as obtained by ND Rietveld rehnement, are reported for shake of comparison. 



Zn K 



Ir Li 


Shell 

N R ^ 

Rnpd 

Shell 

N R ^ 

Rnpd 


(A) (A^) 

(A) 


(A) (A^) 

(A) 

Zn-O 

6 2.05(1) 6E-03 

2.08 

Ir-O 

6 1.96(1) 5.4E-03 

1.9-2.0 




Ir-Ir 

1 2.77(2) 4.1E-03 

2.75 

Zn-Ba 

8 3.55(2) 9.0E-03 

3.58 

Ir-Ba 

7 3.44(2) 1.4E-02 

3.52 

Zn-Ir 

6 3.98(4) l.OE-02 

4.00 

Ir-Zn 

3 4.01(4) 1.2E-02 

4.00 


TABLE 1. (Supplementary Information) Structural data and agreement parameter^ 


Atom 

Site 

X 

y 

Z 

B{K^) 

Occ. 

Bal 

2b 

0 

0 

1/4 

0.23(3) 

1 

Ba2 

4f 

1/3 

2/3 

0.5876(1) 

0.26(2) 

1 

Zn 

4f 

0 

0 

0 

0.28(3) 

1 

Ir 

2a 

1/3 

2/3 

0.15411(5) 

0.34(1) 

1 

01 

6h 

0.4848(2) 

0.9696(5) 

1/4 

0.44(1) 

1 

02 

12k 

0.1701(2) 

0.3402(5) 

0.08318(6) 

0.60(1) 

1 


^ Space group: PGalmmc, a = 5.77034(9) A, c = 14.3441(2) A, V = 413.62(1) A^, = 1.38, Up = 7.04%, = 7.71%, 

Rexp = 6.57%. 



GGk+U 

GGA-bt/-bSOC 


Energy 

(meV/Ir) 

Spin Moment 
(mb) 

Energy 

(meV/Ir) 

Spin (Orb) Moment 
(ms) 



Ir 

01 

02 

Tot 


Ir 

01 

02 

Tot 

FM 

AFMl 

0.0 

44.97 

1.09 

0.78 

0.10 

0.00 

0.19 

0.13 

8.00 

0.00 

0.0 

-41.51 

0.62 (0.10) 
0.73 (0.26) 

0.05 

0.00 

0.12 

0.11 

4.61 

0.00 

AFM2 

-6.59 

0.98 

0.10 

0.14 

0.00 

1.22 

0.65 (0.20) 

0.05 

0.10 

0.00 

AFM3 

26.96 

0.86 

0.00 

0.18 

0.00 

-34.68 

0.74 (0.22) 

0.00 

0.14 

0.00 


TABLE 11. (Supplementary Information) Energy of various magnetic states along with the spin and orbital moments within 
GGA-l-17 and GGA-|-t/-|-SOG approach. Energy for EM state is assumed to be zero. Orbitals moments of Ir are written in 
the parenthesis. 


































